Phonon-induced quadrupolar ordering of the magnetic superconductor TmNi2B2C 
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We present synchrotron x-ray diffraction studies revealing that the lattice of thulium borocarbide 
is distorted below Tq ~ 13.5 K at zero field. Tq increases and the amplitude of the displacements is 
drastically enhanced, by a factor of 10 at 60 kOe, when a magnetic field is applied along [100]. The 
distortion occurs at the same wave vector as the antiferromagnetic ordering induced by the a-axis 
field. A model is presented that accounts for the properties of the quadrupolar phase and explains 
the peculiar behavior of the antiferromagnetic ordering previously observed in this compound. 

PACS numbers: 74.70.Dd, 75.25.+Z, 75.80.+q 



More than ten years ago it was discovered that four 
of the rare-earth borocarbides (Tm, Ho, Er, Dy) shows 
coexistence of superconductivity and antiferromagnetic 
ordering with comparable transition temperatures 1, ^] . 
Superconductivity occurs below = 11 K and the Neel 
temperature is = 1.52 K in TmNi2B2C 3]. The char- 
acterization and the understanding of these intermetallic 
compounds have made much progress during the passed 
decade, see the recent review Although unusual phe- 
nomena have been detected, the type-II superconduc- 
tivity of these materials seems to be described by the 
BCS-theory. The electronic system strongly influences 
the acoustic and optical A4-phonon branches close to 
(0.5, 0, 0) and (0, 0.5, 0) in the non-magnetic Lu and Y 
versions of these compounds 0, S ■ The two vectors 
connect parallel areas of the Fermi surface , so nesting 
explains why the electron-phonon interaction is particu- 
larly large at these wave vectors. This enhancement of 
the electron-phonon interaction is probably the primary 
reason for the relatively high T^-values of the borocar- 
bides. Nesting is generally also assumed to be important 
for the magnetic susceptibility of the band electrons by 
causing a maximum in the Ruderman-Kittel-Kasuya- 
Yoshida (RKKY) coupling of the rare-earth moments at 
the nesting vectors 0. 

The magnetic properties of TmNi2B2C single crystals 
have been studied by magnetization measurements 
and by neutron diffraction 0, Ull H^ . The ordering 
vector at zero field is Qp = (0.094,0.094,0) with the 
magnetic moments transversely polarized along the easy 
c axis. The shift of the magnetic ordering vector away 
from zero is probably a consequence of the Anderson-Suhl 
screening of the long-wavelength susceptibility of the su- 
perconducting electrons Applying a field in excess 
of 10-15 kOe along [100], the magnetic ordering changes 
into another antiferromagnetic one at the wave vector 
Qa ^ (0.483, 0, 0). The new magnetic phase was consid- 



ered to be the stable one at zero field in case the electrons 
were normal, but because the superconducting electrons 
were more strongly affected by the superzone energy gaps 
at than at Qp, the long- wavelength ordering would 
become the most favorable one of the combined system 
at small fields. In the original experiment fll| only fields 
of up to 18 kOe were applied. The extension of the ex- 
periment up to a field of 60 kOe |T3| clearly showed that 
this conjecture could not be correct, since the Q^-phase 
becomes the more stable the higher the field is, and the 
temperature variation more and more resembles that of 
the antiferromagnetic order being induced by the uniform 
field. This highly surprising observation motivated our 
search for a field-dependent quadrupolar ordering of the 
system, observable through the accompanying deforma- 
tion of the lattice. 

The experiments were performed at HASYLAB us- 
ing the BW5 high-energy (« 100 keV) beam-line with 
a triple-axis diffractometer and (1,1,1) Si/Ge-gradient 
crystals as monochromator and analyzer. The single 
crystal of TmNi2B2C, grown by the flux- method 9], is 
plate-like with irregular plate-shape and approximate di- 
mensions 4 X 4 X 0.8 mm'^ along the a, h and c axes, 
respectively. The crystal was mounted with the a and 
c axes in the scattering plane of a horizontal 10 tesla 
cryo-magnet. Geometric constraints required that the 
fleld was offset by 5° from the a axis in order to probe 
the expected lattice distortion wave with the symmetry 
of the A4-phonon mode close to (0.48,0,8) along with 
the (0, 0, 8) Bragg reflection. Temperature and field de- 
pendent superstructure peaks were found at scattering 
vector (0.484(2), 0, 8) with the integrated x-ray intensi- 
ties as shown in Fig. ^ In order to include the intensity 
from minor satellite crystallites we use numerical summa- 
tion. The x-ray cross-section is completely dominated by 
the scattering from the Tm ions, and scans performed at 
(0.484, 0, 0) showed negligible intensity. This is clear evi- 
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FIG. 1: The integrated x-ray scattering intensity at (0.484, 0, 8) as a function of temperature in the presence of an apphed field 
along [100] . The colored symbols used for the different values of the field (in units of kOe) are specified in the legend box. The 
two figures include the same results, but the intensity scales differ by a factor of 10. The intensity is proportional to the square 
of the displacement of the Tm ions, and the solid curves show the results derived from the model. The insert shows the field 
dependence of the strain amplitude £13 at 1.7 K derived from the experiments in comparison with the calculated results. 



dence for a transverse displacement of the Tm-ions along 
the c axis corresponding to the acoustic A^-phonon sym- 
metry. The (0, 0, 8) Bragg reflection was used as reference 
for accurate determination of the superstructure peak po- 
sition and estimates of the superstructure displacement 
amplitudes. However, the irregular crystal geometry and 
the small scattering angles result in rather low and un- 
certain sample transmission coefficients of the order of 
4.6 X 10-3 and 8.0 x lO^^ for the (0, 0, 8) and (0.484, 0, 8) 
reflections, respectively. Including these corrections and 
defining the strain of the Tm ion at site 

Ei,i^) = = i?i3 cos(Q^ ■ fi. + 0) , (1) 

a/ 2 

we find that the absolute displacements of the Tm ions 
at 1.7 K and 60 kOe, correspond to a strain amplitude 
« 6 X 10~3. The resolution of the diffractometer 
was determined by the peak profiles of the (0, 0, 8) Bragg 
peak yielding a crystal mosaicity of w 0.02° FWHM. The 
scans along the a axis were essentially rocking curves 
with Lorentzian profiles while scans along c were Gaus- 
sian shaped. At low temperatures the correlation length 
along a, corrected for instrumental resolution, is deter- 
mined to be 185 A when the quadrupolar phase co-exists 
with the magnetic Qp phase, and it attains an almost 
constant value of ~ 275 A in the phase. At zero 
field and low temperatures the Gaussian peaks in the c 
direction indicate box-like structures of lengths w 230 A 
that increases to values above the resolution limit of ~ 
300 A at 100 kOe. The superstructure peak widths are 



essentially temperature independent at low temperatures 
but start to broaden when approaching Tq. 

Setting up a model we use the crystal-field (CF) pa- 
rameters of the Tm ions determined from the observation 
of CF transitions and from the high-temperature sus- 
ceptibilities 1^ . Here we take into account the anisotropy 
due to the classical dipole-dipole interaction, ^^^(O) = 
3cc(0) + 9-15 /^eV, adjust the CF-parameters accord- 
ingly {B^ = -0.1, Bl = 0.00033, Bt = -0.01, = 
0.78 X 10-^ and S| = -1.1 x lO"'' in units of meV) and 
include the demagnetization field in the Zeeman term. 
The zero-field value of T-^ determines the effective two- 
ion coupfing to be ScdQp) — ^-"^ M^V, which is also 
approximately its value at zero wave vector in the nor- 
mal phase. The purely magnetic part of the Hamilto- 
nian is well characterized, and the additional quadrupo- 
lar/magnetoelastic part of the mean-field Hamiltonian is, 
by symmetry, considered to be 

A^i) - - B,3E,,{i)Ol{i) + CpElS) + iBpEU^ 
-K{QA){0l{i))[0l{i)-\{0\{^)]. (2) 

Here 0\{i) — ^{J^Jz + JzJx)i ^itl^ x, y, and z along 
[100], [010], and [001] respectively. E-^^{i) is determined 
by minimizing the strain-dependent part of the average 
free energy: 

AF = -B,.,E,,{{Ol)) + \ceEI^ + ^BpEt, , (3) 

where ((O2)) is the thermal value (02(«)) cos(Q^-ilj-)-0) 
averaged over the lattice. E^^ determines E^^{i) accord- 
ing to eq. The transition temperature Tq of the 
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quadrupolar ordered phase is independent of the value of 
the fourth-order elastic constant B^. It is determined by 
the effective quadrupolar coupling 

K^siQA) = K{Qa) + ^ (4) 

being equal to the inverse of the non-interacting 
quadrupolar susceptibility at the temperature Tq, or 
K^f^{QA) = 1/xqi (Tq) = 0.0187 meV, when assuming 
Tq = 13.5 K at zero field. The properties of the CF 
Hamiltonian imply that the quadrupolar susceptibility 
increases with the a-axis field, and the predicted field 
dependence of Tq is in reasonable agreement with obser- 
vations according to Fig.^ In the zero-temperature limit 
Xqi becomes even more sensitive to the field. Here its 
maximum value at a field of 20 kOe is nearly a factor of 
3 larger than the zero-field value. Although this increase 
is large, it is far from explaining the dramatic enhance- 
ment of E^^ observed at 1.7 K, see the insert of Fig. ^ 
In the non-magnetic borocarbides, the elect ron-phonon 
interaction is close to enforce a complete softening of the 
A4-phonon mode at the nesting wave vector correspond- 
ing to Q A in the Tm system. The phonon energy is 10 
meV at room temperature and about 4 meV at 4.2 K 
in the Lu compound 1^. In the model above, the ef- 
fects of the electron-phonon interaction are included in 
terms of the phenomenological parameters and B^. 
If the electron-phonon interaction is as important in the 
present system as in Lu borocarbidc, then should be 
small, nearly zero, and B^ would be decisive as soon as 
the lattice distortion becomes non-zero. This is the argu- 
ment for introducing _B^, which has the additional effect 
that the two contributions to the effective quadrupolar 
coupling behave differently. This allows a differentia- 
tion between the effects due to B^^^ and those due to 
the quadrupolar-quadrupolar interaction K{Qa)- 

The quadrupolar ordering has the unique consequence 
(we are not aware of any other systems showing this be- 
havior) that the application of a uniform field induces a 
modulated magnetic moment perpendicular to the field. 
This effect was the very reason why we started to look 
for the quadrupolar ordering, and the explanation for 
it is straightforward. When £'13(1) or {0\{i)) is non- 
zero, the CF Hamiltonian includes an anisotropy term 
0\{i). This implies that the paramagnetic moment of 
the ith site induced by an cc-axis field is not along the 
field, but is rotated an angle towards the z axis in pro- 
portion to the size of the anisotropy term. Since the 
anisotropy, proportional to {0\(i)), is modulated in a 
sinusoidally way with the wave vector Q^, the z compo- 
nent of the induced moment is going to show the same 
modulation. The tetragonal symmetry of the system 
implies that the Hamiltonian also includes the equiva- 
lent terms obtained by interchanging x and y in eq. l(2Jl. 
In the domains in which {02^{i)) is ordered with the 



wave vector = (0,0.484,0), the application of a field 
along [100] does not produce any perpendicular magnetic 
component, and the quadrupolar order parameter is only 
weakly changed by the field. This means that the [100] 
field induces a modulated z-axis moment at the wave 
vector Qj^ but not at Q^, and it destabilizes the 
domains. Hence, the first circumstance explains the ab- 
sence of mag netic scattering peaks at iQ^ discussed in 
Refs. [lllll2j . In the final comparison between the model 
and the experiments, we have introduced one more pa- 
rameter, namely, the value of '5cc{Qa)- This parameter is 
not decisive for the magnitude of {Jc{Qa)) (notice that 
Tq ^ TJsr)' t)ut is important for E^.^ at the lowest tem- 
peratures, and the fit is significantly improved, when as- 
suming this parameter to be small or, actually, slightly 
negative. The final model is 

IB13I = 1.91, ce = 360, Be = 4.8 x 10^ 
K{Qa) = 0.0086, dcAQA) = -0.0029 (5) 

all in units of meV. The different quantities have been 
calculated, see Ref. Hi, by assuming a commensurable 
quadrupolar ordering with a period of 25 lattice planes 
along the a axis, i.e. Qa = (0.48, 0, 0). The square of the 
distortions is observed to be a factor of 2-3 smaller than 
calculated at fields below 15 kOe. An equal population 
of the Qj^- and Q^-domains, at the lowest values of the 
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FIG. 2: The temperature dependence of the integrated neu- 
tron scattering intensity at Q 4 = (0.484, . 0) at various val- 
ues of the applied field along [100], see Ref. 1121 . The different 
values of the field (in units of kOe) are specified in the box. 
The results were obtained in two different experimental set 
ups (circles and squares). The repetition of the experiment 
at 40 kOe was used for determining a common scale for the 
two sets. The intensity is proportional to the squared ampli- 
tude of the transverse magnetic moment at the wave vector 
Qj^ and the solid lines are the calculated results. 
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field, may account for a factor of 2. The model neglects 
any direct consequences of superconductivity on the mag- 
netic phase diagram, but the experimental x-ray scatter- 
ing intensities show kinks in their temperature variations 
at the points where the field is equal to the supercon- 
ducting critical field H^2- This indicates a (moderate) 
reduction of the quadrupolar order parameter in the su- 
perconducting phase. 

The neutron-diffraction results [l^ are compared with 
theory in Fig. El assuming that the value 6 on the inten- 
sity scale corresponds to the moment amplitude of 5.0 
/i^ . The calculated amplitudes of the modulation of the 
moments are of the right order of magnitude. The max- 
imum of the amplitude in the zero-temperature limit is 
predicted to occur as soon as the system leaves the Qp- 
phase, at a calculated field of 17 kOe to be compared 
with an experimental value of 10-15 kOe. The experi- 
mental maximum is at about 40 kOe, but wc have to 
remark that there are relatively large, unexplained scat- 
ter in the experimental results. Instead of focussing on 
the discrepancies we rather emphasize the ability of the 
present model to account for the highly disparate varia- 
tion of the properties of the system. Using a simple elas- 
tic model we estimate that = 360 meV corresponds 
to an acoustic phonon mode at 2.4 meV close to the soft 
phonon energy of 4 meV observed in the Lu compound. 
The value of also appears to be acceptable. The en- 
tropy of the electron-phonon system has nearly vanished 
and the electron-phonon parameters and are ex- 
pected, as assumed, to stay practically constant below 20 
K. Because of the small non- magnetic entropy, the heat 
capacity anomaly at Tq is expected to be weak, since 
the CF contribution of the 4/ electrons is negligible (the 
ground-state doublet remains degenerate at zero field). 

The analysis indicates that the electron-phonon sys- 
tem in TmNi2B2C, like in the Lu compound, is very close 
to the point of a lattice instability. In TmNi2B2C, the 
transition is occurring because of the additional contribu- 
tions of magnetoelastic and quadrupolar-quadrupolar in- 
teractions. The magnetoelastic interaction derives from 
the electrical crystal field produced by the displacements 
of the ions, whereas the quadrupolar-quadrupolar inter- 
action may be mediated by the same charge-density wave 
responsible for the strong electron-phonon interaction, 
and, hence, may be an additional consequence of the 
Fermi-surface nesting. One of the surprises of the present 
analysis is that the RKKY-exchange interaction does not 
seem to be much influenced by the nesting, since the in- 
teraction is found to be much smaller at the nesting wave 
vector than at Qp. 

In conclusion, the x-ray experiments exposed the pres- 
ence of a quadrupolar phase in TmNi2B2C below 13.5 K 
at zero field. The transition temperature increases and 
the accompanying deformation of the lattice is grossly 
enhanced when a field is applied along [100]. From the 
model analysis we conclude that this occurs mainly be- 



cause the lattice is already close to the critical point due 
to the non-magnetic electron-phonon interaction. The 
magnetic interactions are the final ingredients that en- 
force the system to enter the quadrupolar phase, and 
the analysis indicates that the magnetoelastic crystal- 
field and the quadrupolar-quadrupolar interactions are of 
nearly equal importance, K{Q^) = 0A6K^g{Qj^). The 
quadrupolar ordering has the unique effect that a field 
applied along [100] induces antiferromagnetic ordering at 
the same wave vector as the quadrupolar ordering. The 
Fermi-surface nesting is presumably important for the 
strong electron-phonon interaction at Q^, and possibly 
also for the quadrupolar-quadrupolar interaction. Con- 
trary to this, the RKKY interaction is concluded to be 
relatively weak at this wave vector. This suggests that 
those electrons, which are close to the nesting areas and 
important for creating the superconducting phase, are 
not the same as those determining the RKKY interac- 
tion. The present discovery has basic consequences for 
the understanding of the magnetism of the borocarbides. 
It raises a number of questions: Is there a strain/charge 
quadrupolar phase in the non-magnetic borocarbides or 
a magnetic quadrupolar phase in other of the magnetic 
ones? Is the strong peak at (0.55, 0, 0) in the exchange 
interaction of the Er compound, see Refs. lilEl, a con- 
sequence of a quadrupolar interaction rather than the 
RKKY? Is the quadrupolar interaction the reason for the 
longitudinal polarization of the antiferromagnetic or- 
dered moments in the Tb compound? etc. 
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